#include "cHydrogen.h"

Hydrogen::Hydrogen()
{
	setName("Hydrogen");
	setUniverse(30.0);
	setWidth(getUniverse());
	setPotLeft(0.0);
	setPotRight(getUniverse());
	setFx(getPotRight());
        setUniLeft(0.0);
        setUniRight(getUniverse());
	setTx(0.0);
	setMinEnergy(-1.0);
	set_l(0);
	//filter[activefilter]->calculate(this);
}

Hydrogen::~Hydrogen()
{
}

double Hydrogen::getPotentialValueAt(double x) {
    if(x<=0.0001) {
        return -10000.0;
    }
    else {
	return -1.0/x + (get_l()*(get_l()+1)/(2*x*x));
    }
}
